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Calculations of the first and second 
moments of displacement damage 
energy distributions from clastic colli- 
sions and from nuclear reactions, at 
proton energies ranging from 10 MeV 
to 300 MeV, are incorporated into a 
model describing ihe probability of 
damage as a function of the proton 
fluence and the size of the sensitive 
micro-volume in Si. Comparisons be- 
tween the predicted and measured 
leakage currents in Si imaging arrays 
illustrate how the Poisson distribution 
of higher energy nuclear reaction 
recoils affects the pixcl-to-pixel variance 
in the damage across the array for pro- 
ton exposures equivalent to mission 
duration of a few years within the 
earth's trapped proton belts. Extreme 
value statistics (EVS) quantify the 
largest expected damage extremes 



following a given proton fluence, and 
an analysis derived from the first-princi- 
ple damage calculations shows excellent 
agreement with the measured extremes. 
EVS is also used to demonstrate the 
presence of high dark current pixels, or 
"spikes," which occur from different 
mechanisms. Different sources of spikes 
were seen in two different imager 
designs. 
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1. Introduction 

Proton-induced displacement damage degrades 
semiconductor electrical properties by introducing 
localized energy states within the band-gap which 
result in increased generation dark current, carrier 
recombination and charge trapping. On average, 
the permanent proton-induced damage in bulk Si is 
proportional to the average amount of energy 
which has been imparted through non-ionizing 
processes following elastic and inelastic scattering 
of Si atoms [1-3]. However, on micro-volume scales 
appropriate for microelectronics, average damage 
is a poor indicator of damage effects because of 
differences in the number of incident particles and 



fluctuations in energy deposition which are an 
unavoidable consequence of collision kinematics. 
Characterization of displacement damage in Si 
micro-volumes has particular importance for satel- 
lite imaging array applications. Device radiation 
hardening solutions have largely solved problems 
associated with ionization effects. However, parti- 
cle irradiation seriously degrades charge transfer 
efficiency through carrier trapping and increases 
dark current by carrier generation. Permanent dark 
current increases from single particle interactions 
have been reported in sensor arrays following pro- 
ton and neutron irradiation [4,5]. Pixel-to-pixel 
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variations in dark current increases following 
multiple interactions within each pixel have also 
been shown to depend on the incident particle and 
energy [3,6]. 

Orbital proton energy spectra, whether from the 
earth's trapped radiation belts or solar flares, typi- 
cally peak at very low (and more damaging) ener- 
gies and decrease exponentially with increasing 
proton energy. Typical spacecraft structural shield- 
ing effectively attenuates lower energy protons re- 
sulting in spectra extending from a few MeV to 
several hundred MeV with average energies over 
20 MeV. Proton linear accelerators and cyclotrons 
are therefore well suited for monoencrgetic charac- 
terizations of damage verses proton energy which 
can then be incorporated in damage predictions for 
a given environment and shielding configuration. 

For the proton energy range of 10 MeV to 
300 MeV, this work explains the average damage 
and pixel-to-pixel damage fluctuations in terms of 
calculated parameters reflecting the energy depen- 
dence of the proton-silicon interactions. The analy- 
sis predicts the damage distribution within a given 
array as illustrated for the particular case of a 
charge injection device (C1D) depletion volume 
and the cross-sections and Si recoil energies 
applicable to 12 MeV, 22 MeV, and 63 MeV pro- 
ton induced damage. This enables a direct com- 
parison between the predicted damage distri- 
butions and the observed dark current histograms 
reported in [3] across a range of energies important 
for orbital environments. 

Dark current extremes, which may follow from 
damage extremes, are a particularly serious con- 
cern for a variety of satellite imager applications. 
These "hot pixels" of "spikes" interfere with the 
instruments ability to resolve small, dim objects 
such as low magnitude stars which might be used 
for a star tracker guidance system. Also, spikes in a 
image can place overhead on data compression al- 
gorithms and burden telemetry channels. Extreme 
value statistics are well suited for characterizing 
the frequency and magnitude of these spikes. 
These tools are applied to proton damaged CID 
imagers to illustrate this approach, and we show 
that for one particular CID design, the spikes can 
be accurately predicted based on the calculated 
probabilities and kinematics of proton-initiated 
nuclear reactions. 



2. Recoil Spectra Parameters 

For proton energies of practical interest in satel- 
lite orbits, the damage is caused by recoiling atoms 
from collisions with Si atoms. As depicted in Fig. 1 
[7], elastic scattering by the Coulombic field of the 
nucleus dominates for protons below 10 MeV, 
though at higher energies, nuclear elastic scattering 
also becomes important. By 60 MeV, about half of 
the displacement damage is due to nuclear inelastic 
reactions which dominate above 100 MeV. Elastic 
cross-sections are relatively high with recoil ener- 
gies typically less than 1 keV as opposed to in- 
frequent nuclear reactions emitting very damag- 
ing MeV-range recoils. In this work, the first and 
second moments of the recoil spectra are calcu- 
lated separately for each type of interaction. 

The average damage energy from all elastic 
recoils is obtained by numerically integrating the 
product of the differential cross section weighted 
by the corresponding recoil damage energy, over all 
scattering angles. Damage energy is defined here 
as the portion of energy lost by a recoil through 
mechanisms other than ionization as calculated by 
Lindhard et al. [8]. Note that this represents an 
important adjustment to the total energy imparted 
by the reaction atoms which must be assessed for 
evaluating either the nonionizing or the ionizing 
energy imparted. The second moment calculations 
proceed in the same manner, except the recoil 
damage energies now appear to the second power. 
The variance follows as the second moment minus 
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Fig. I. Frequent coulombic scattering from protons Of a 
few MeV initiate low energy recoil atoms resulting in isolated 
defect sites. More energetic protons can impart more energy to 
recoil atoms via nuclear elastic and inelastic reactions resulting 
in less frequent but mure complex damage structures. 
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the first moment squared, as is customary. Figure 2 
plots the mean and variance of non-ionizing energy 
for proton energies from 10 MeV to 300 MeV, and 
Table 1 lists the values of experimental interest, 
along with the total elastic cross-sections and recoil 
energies. 

The inelastic reaction cross-sections are esti- 
mated according to the empirical formula of Letaw 
et el, [9]. Calculations of primary recoil energies 
consider both the initial intranuclear cascade and 
subsequent evaporation of nucleons. The momen- 
tum imparted during the evaporation phase is esti- 
mated using a Brownian motion model. Next, the 
average and variance of the damage energy are 
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Fig. 2. The mean and variance of ihc lotal elastic damage energy 
arc plotted versus proton energy along with a bcsi-fit curve. The 
moments were calculated based on clastic differential cross- 
section data [1 and references therein] indicated by circles and 
triangles. 



calculated as in the elastic case, and the results are 
summarized in Table 1. Further details and 
comparisons with data are discussed in [1,3]. 



3. Damage Calculations 

For a given proton energy, the mean and variance 
describing the probability density function (pdf) 
for damage from single interactions, as listed in 
Table 1, allow independent evaluation of the dam- 
age expected from the elastic and inelastic recoil 
categories. For the elastic category the mean for the 
pdf describing damage at a given proton fluence is 
the product of the number of interactions and the 
mean of the pdf for single interactions. The number 
of interactions is the product of the average cross- 
section, the incident particle fluence, and the num- 
ber of Si atoms in the sensitive volume. 

The elastic scattering component of the variance 
associated with the fluence dependent pdf is esti- 
mated as the product of the number of interactions 
and the single interaction pdf variance shown in 
Fig. 2. This is possible because Poisson fluctuations 
in the number of elastic recoils per pixel do not con- 
tribute significantly to the final result. In the regime 
where N, the average number of interactions per 
volume element, is greater than 20, the Af-fold con- 
volution of the single interaction pdf with itself 
leads to a Gaussian elasiic damage distribution with 
mean and variance as described above. 

For sensitive volumes and fluences of interest 
here, the average number of inelastic recoils ranges 
typically from a fraction to a few, and a discrete 
Poisson distribution determines the probability of a 
given number of inelastic recoils. The pdf governing 
the inelastic damage energy for a pixel with yVj 



Tabl« I. Proton recoil spectrum parameters 
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22 
63 
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0.569 
1.44 
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0.152 


2.05 x t0 - 3 
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inelastic recoils reflects the /Vj-fold convolution of 
the pdf for single inelastic damage. For purposes of 
this analysis, the form of the single event pdf for 
inelastic recoil products is approximated as a two 
parameter gamma distribution with mean and vari- 
ance as indicated in Table 1. 

Since the elastic and inelastic processes are in- 
dependent random variables, the combined damage 
for pixels in which both occur follows as the con- 
volution of the pdfs describing each of the two 
categories. Figure 3 illustrates this simulation for 
the specific case of the imaging array used in this 
study in which damage from a fluence of 4.0 x 10 H1 
12 MeV protons/cm 2 occurs, and each pixel's sensi- 
tive volume is 1300 jim 3 . The Gaussian distribution, 
shown in Fig. 3a as the N, =0 case, describes dam- 
age corresponding to an average of 4,000 events per 
pixel. Figure 3a curves for A/} = 1 through 10 inelas- 
tic recoils per pixel reflect increases in both the 
means and variances as the shape tends toward 
Gaussian. Figure 3b shows the pdf for total com- 
bined damage as the superposition of the pdfs in 
Fig. 3a, after weighting by their associated Poisson 
probabilities according to the average of 1.8 inelas- 
tic recoils per pixel. This average is arrived at by 
considering the number of silicon atoms present in 
the 1300 u.m 3 volume, and the composite cross- 
section for nuclear inelastic reactions for 12 MeV 
protons as shown in Table 1. 
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Fig. 3 (a). The Gaussian distribution with no inelastic recoils de- 
scribes elastic damage and convolved distributions show com- 
bined damage from clastic and 1 through 10 inelastic recoils, 
(b) Weighting according tn the Poisson probabilities precedes 
the superposition to determine combined damage probabilities. 
The simulation applies to the Si CTD sensitive volume of 
1300 u.m v and 4.0x 10'" 12 MeV protons/cm 1 . 



Early in a space mission or in a relatively benign 
orbit, the fluences may be 1-2 orders of magnitude 
lower, at about 10 s cm -2 . In this low fluence 
regime, the very low probability of inelastic recoils 
suggests that two would probably not be observed 
in the same volume element. The number of elastic 
recoils per volume would be correspondingly low 
resulting in very large relative changes within the 
pixels where nuclear reactions occur. The product 
of the low probability of an inelastic event with the 
large number of pixels determines the pixel popu- 
lation for which damage exceeds the average by 
factors of up to 1,000. 



4. Predicting Damage Extremes 

In addition to being a necessary tool for assessing 
radiation-induced fixed pattern noise, the probabil- 
ity density function describing damage throughout 
the array can be used to predict the number of 
elements sustaining exceedingly large damage in- 
creases afler a specified exposure. In [6] it was 
shown how extreme value statistical analysis can 
describe the measured distribution of pixels with 
the largest damage increases following 12 MeV and 
63 MeV proton damage to the Si CID. For a broad 
range of proton energies and fluence levels, the 
largest extremes were shown to obey a Type 1 ex- 
treme value distribution. Next it will be shown that 
the particular Type 1 distribution describing pro- 
ton-induced damage extremes can be predicted 
from the calculated pdf described above. 

Figure 4 shows an expanded view of the tail 
region in Fig. 3b which identifies the contributions 
to the pdf from the 11 populations containing 
through 10 inelastic recoil events per pixel. The 
damage energy distribution has a mean of 
0.85 MeV, and the skewed high energy tail extends 
to about 1.8 MeV. Individual distributions are iden- 
tified according to the number of inelastic recoils. 
Figure 4 illustrates how several of the component 
distributions contribute to the probability of ex- 
ceeding large damage energies. Based on a total 
pixel population of 61,504, the inset presents the 
number of pixels expected above the specified 
damage level, E d . This number is the total popula- 
tion multiplied by pj, the probability of exceeding 
damage energy E& within the whole array. This 
probability is calculated as the summed pdf inte- 
grated from Et to infinity. 
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Fig. 4. For the simulation depicted in Fig. 3, volumes containing 
from 1 to 10 inclasiic recoils contribuitc to the population of 
pixels with the most damage. The inset show's the number 
expecied above a given damage energy for a 61,504 pixel array. 



Two steps are necessary to compare these results 
on the basis of the cumulative Type 1 extreme value 
distribution. As discussed in [6,10-12], extreme 
value analysis can be applied to data to evaluate the 
probability of exceeding a certain value within any 
population size by evaluating a set of largest values 
extracted from subsets of a given population. In the 
next section we will treat the case where the 61,504 
pixel population has been subdivided into 248 
groups of 248 pixels each. Using p d as defined 
above, the probability, po, of having no pixels ex- 
ceeding Ei within the group of 248 pixels can be 
evaluated using the discrete binomial distribution 
as: 



po - (!-/><)) 



248 



(1) 



The associated standard variate specific to the ex- 
treme value cumulative probability plot is given by: 



S(p») - -ln[-ln(p )]. 



(2) 



Thus £<j, or a proportional quantity such as dark 
current, can be plotted against the corresponding 
standard variate to predict the Type 1 extreme dis- 
tribution specific to the pdf from which it is gener- 
ated. Detailed discussions of extreme value analysis 
are discussed in the references [11,12], and applica- 
tions to this study will be illustrated in the following 
section. 

5. Comparison with Dark Current Data 

Calculations described in the previous section are 
compared here to measured dark current increase 
distributions specific to proton-induced damage in 
a General Electric 256 pixel x 256 pixel Si CID. 



Devices are fabricated in an n-type Si epitaxial 
layer doped with 5 x 10 14 P atoms/cm 2 . A field iso- 
lation oxide confines the collection area to about 
17 mm x 17 mm, but for purposes of dark current 
studies only the 1300 u.m 3 depletion volume leads 
to carrier generation. 

All dark current data reported here were ac- 
quired at 18.0 °C and correspond to a 248 x 248 sub- 
set of the array. After each proton exposure and 
measurement the dark current increase for each 
element was calculated by a pixehby-pixel subtrac- 
tion of the pre-irradiation value. This correlation 
removes imager spatial noise not resulting from 
radiation. Temporal read-out noise accounts for 
less than 5% of the dark current spreads reported 
here. More detailed descriptions of the imaging 
array and the dark current measurement are pro- 
vided in [6]. 

Proton irradiations with energies of 12 MeV, 
22 MeV, and 63 MeV were performed at the 
University of California at Davis cyclotron facility. 
The beam line and dosimetry have been described 
previously [13]. Irradiations were conducted at a 
nominal dose rate of 1 kRad(Si)/s with all leads 
grounded. Dark current measuremenls were ini- 
tiated about 15 minutes post irradiation and re- 
peated after I day and again after about 1 week. No 
significant annealing was observed over this period. 

In Fig. 5, comparisons are made between dark 
current data histograms and calculated damage 
energy distributions in the CID pixels. The calcula- 
tion approach described above has been exercised 
for three 12 MeV proton fluences corresponding to 
averages of (1.8, 4.5, and 9.0) inelastic recoils/pixel. 
Based on the population of 61,504 pixels and 
Poisson statistics, the maximum numbers of inelas- 
tic recoils expected in any single pixel are 10, 16, 
and 24, respectively. For comparing the calculations 
to dark current data, a conversion factor relates the 
average dark current and the mean damage energy. 
For the three fluences, the average conversion fac- 
tor of 2.2 nA/cnr per MeV of damage energy varies 
by up to 10%, which reflects the experiment's 
dosimetry uncertainty. The calculated damage 
curves in Fig. 5, based on the first and second 
moments of the non-ionizing energy imparted by 
the recoil spectrum, describe the dark current dis- 
tribution to a remarkable degree of accuracy. 

Comparisons for 22 MeV and 63 MeV proton 
damage show similar agreement. The coefficient of 
variation, defined as the ratio of the standard 
deviation to the mean damage (or dark current), is 
a dimensionless figure-of-merit. At 12 MeV, 
22 MeV, and 63 MeV, the experimental and calcu- 
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Fig. 5. Calculated damage energy distributions show excellent agreement with 
measured dark current histograms from a Si CID damaged by 12 MeV proteins. 
Calculations arc based on averages of (1.8, 4 J, and 9.0) Inelastic recoils per pixel, 
and the damage distribution shapes reflect the associated discrete Poisson distribu- 
tions. 



lated results agreed within 2%, 12%. and 15% 
respectively [14]. Also at 63 MeV, wilh 45% of the 
damage caused by inelastic recoils, the means of 
the two distributions are normalized by a factor of 
2.0 nA/cm 2 per MeV of damage. This does not 
differ significantly from the conversion factors de- 
termined for 12 MeV, thus demonstrating that the 
average damage is proportional to the energy lost 
through non-ionizing processes, and that the ex- 
pected damage from both the elastic and inelastic 
categories is present. 

The somewhat better agreement between calcu- 
lated and measured damage distributions at the 
lower proton energy of 12 MeV could be influenced 
by characteristics associated with high energy 
recoils. At proton energies of 12 MeV and 63 MeV, 
the contribution to the total damage from inelas- 
tic reaction recoils increases from roughly 15% to 
45%. Also, as this fraction increases, the average 
inelastic recoil energy (and range) also increases, 
and at higher proton energies the higher energy 
recoil ranges approach the smallest dimension 
of the sensitive volume (about 2 u.m). These issues 
would be even more important for smaller sensi- 
tive volumes (i.e., CTE loss in a CCDs buried 
channel). 



6. Largest Dark Current Extremes 

Here the measured largest dark current in- 
creases are compared to the calculated damage 
maxima for the specific cases of the three 12 MeV 
proton fluences of Fig. 5. For each proton energy 
and fluence level, the dark current extreme popula- 
tions are generated by subdividing the 61,504 pixel 
population into 248 groups of 248 pixels each. The 
largest value from each group comprises the popu- 
lation of extremes. Figure 6 depicts how the ex- 
treme distribution is derived for the case of the 
lowest fluence level shown in Fig. 5 (note this 
example also corresponds to the calculations for 
Figs. 3b and 4). After ranking and estimating the 
probability according to the [rank/(n+l)] for n 
samples as in [6], the standard variate follows from 
Eq. (2), and the measured dark current extremes 
can be compared with the Type 1 extreme value 
distribution using a Type 1 cumulative probability 
chart. 

Likewise, damage maxima calculated as de- 
scribed in section C can be compared to the same 
Type 1 extreme probability distribution using 
Eqs. (1 and 2) and the normalization constant of 
2.2 nA/cm 1 per MeV of damage energy. Figure 7 
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compares measured dark current extremes, for the 
three 12 MeV proton fluences treated in Fig. 5, to 
predicted damage maxima according to the Type 1 
distribution. The linear character of the data and 
calculation show that they obey a Type 1 distribu- 
tion, and the close agreement at each fluence 
demonstrates the robustness of the analysis. The 
return period abscissa at the top of Fig. 7 identifies 
the largest expected dark current increase for a 
given number of array subsets. For example, at the 
fluence of 2.0 x 10'Vcm 2 the return period value of 
10 corresponds to about 13 nA/cm : indicating the 
largest expected increase within a set of 10 groups 
or 2,480 pixels. Good agreement also exists be- 
tween the measured and predicted extremes from 
63 MeV protons. 



The ability of the calculation to predict the 
largest measured dark current changes offers in- 
sight into the mechanisms responsible for proton- 
induced damage extremes. The linear response on 
the Type 1 plot indicates that a single mechanism is 
probably responsible for largest values while the 
slope reflects the variance. As pointed out in Fig. 4, 
the largest damage regions in this fluence regime 
follow from the probabilistic treatment of pixel 
populations damaged by several inelastic reaction 
recoils. 

When the probability of an inelastic recoil per 
pixel is much less than one, as is the case in many 
natural space environments, the analysis can deter- 
mine the total number of pixels expected with dam- 
age above a given level. In this regime, where the 
background radiation-induced damage can be quite 
low, largest damage regions can be several hundred 
times the average. Some of the array subsets would 
have largest changes dominated by single inelastic 
recoil damage and others by the largest of the less 
damaging elastic recoils. In this case, agreement 
with the Type I cumulative chart could be expected 
only with sufficiently larger bin sizes so that each 
bin would include at least 1 pixel with damage from 
an inelastic reaction. 

A qualitative comparison of such a situation fol- 
lows from our evaluation of the proton response of 
an alternate CID imager design. The important as- 
pects of this "narrow row" design were previously 
discussed in [6], with the key difference resulting in 
spurious high electric field profiles near the row 
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Fig. 7. Cumulative probability distributions demonstrate excellent agreement between 
calculated damage extremes and the measured dark current extremes b3scd on a 24S 
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electrodes. The comparison of the dark current 
and extreme distributions for this device type, 
shown in Fig. 8, can be made with the previously 
discussed design at the same proton exposure level, 
as in Fig. 6. Note that the average dark current is 
doubled, but more importantly, the character of 
the extreme distribution is markedly different. The 
consequence of this is more evident in the proba- 
bility chart of Fig. 9. Clearly the narrow row design 
results in an extreme distribution which is not Type 
1 when analyzed as before. Rearrangement of the 
array to 31 bins of 1984 pixels offered a better 
match with the Type 1 distribution. Even so, the 
extremes for this case cannot be understood based 
on first principles analysis of damage mechanisms 
as before. We concluded that in this case, the 



largest extremes were not caused by conventional 
charge generation, and extreme value statistics 
played a critical role toward quantifying the likeli- 
hood and magnitudes of this other mechanism. In 
[4] we discuss supplemental measurements and 
analysis which have lead us to conclude that the 
high field regions were causing localized lowering 
of the band-gap resulting in field enhanced emis- 
sion and tunneling currents. Thus the statistics of 
extremes are applied to evaluate design variations 
and to assure that optimum imager performance 
can be assured. We also concluded that acceptable 
designs should have extreme characteristics as 
depicted in Fig. 6 which are limited only by un- 
avoidable consequences of particle-semiconductor 
physics. 
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suggests the role of field enhanced mechanism in causing the largest leakage 
extremes. 



7. Conclusions 



This paper presents an analytic approach for de- 
termining the pixel-to-pixel distribution of particle- 
induced displacement damage in micro-volumes 
representing sensitive volumes in sensor arrays. 
The calculation is based on interaction cross- 
sections as well as parameters describing the dam- 
age imparted by the spectrum of particle-initiated 
recoils. It predicts the dark current distribution 
and largest dark current changes in a Si CID fol- 
lowing incremental damage with 12 MeV, 22 MeV, 
and 63 MeV protons. These proton energies span a 
regime important to the natural space environ- 
ment; lower energy protons for which the damage 
is dominated by elastic scattering and higher ener- 
gies where nuclear reactions become increasingly 
important. The analysis illustrates how high energy 
recoils from nuclear reactions influence the pixel- 
to-pixel variance in proton-induced damage and 
cause the largest damage occurrences. To under- 
stand the important exception, we rely on extreme 
value statistics to identify and quantify the role of 
electric field enhanced emission as a mechanism 
for excessive leakage currents. 

The calculation is general in the sense that once 
the parameters describing the recoil spectrum are 



determined, the particle-induced damage distribu- 
tion can be calculated as a function of particle type, 
particle fluence, sensitive volume, and material. 
The significance of these results is that once the 
factor relating the mean dark current to the dam- 
age energy is known from a single measurement on 
a particular array, the radiation response in a 
specified environment can be predicted. In addi- 
tion to providing a means for assessing the radia- 
tion response of a given imager, this analysis has 
flexibility enabling the design-phase evaluation of 
the radiation response of different pixel geometry 
and materials in a variety of environments. 

Extreme value statistics play a critically impor- 
tant role in understanding leakage current spikes 
and in assuring reliable satellite performance. In 
ongoing related research we continue to rely on 
this valuable tool for assessing damage and single 
particle ionization extremes in infrared imaging 
arrays and in optoelectronic detector materials for 
high data rate spacecraft data links, each of which 
must perform to exacting standards to assure reli- 
able performance of extremely valuable space 
assets. 
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